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1. Introduction 

In the past two decades the study of chaotic dynamical systems has attracted 
attention of many physicists. Roughly speaking, a system is said to be chaotic 
if orbits of the motion are in some sense irregularly distributed. One of the 
most useful measures of this irregularity is the Kolmogorov-Sinai entropy. 
Another important quantity relevant to control the chaotic behavior is the 
Lyapunov exponents, which measure the exponential instability of almost all 
orbits with respect to the change of initial conditions. It turns this instability 
leads to the loss of memory of initial conditions, decay of correlations and 
approach to statistical equilibrium. 

It is believed that also quantum systems have different qualitative proper- 
ties depending on whether the corresponding classical systems are integrable 
or chaotic [1]. For example, they reveal significant differences in the char- 
acter of their wave functions or the distribution of their energy levels [2]. 
Also the numerical analysis of some finite quantum systems shows quantita- 
tive differences in regular and chaotic regime [3]. It is obvious that on the 
quantum level the distinction between chaotic and quasiperiodic behavior 
must be blurred, nevertheless some relative measure of the degree of chaos 
should exist. There have been many attempts to find fingerprints of chaotic 
behavior in quantum dynamical systems. The most natural one is based on 
the correspondence principle, which states that a quantum system is chaotic 
if its classical limit is chaotic. However, it can only indicate some features 
of chaotic behavior but cannot serve as a precise definition. As was pointed 
out by van Kampen [4] such notions like ergodicity or mixing need a limit 
t — > oo while the correspondence principle refers to /i — > 0, and these two 
limits do not commute. Moreover, an example of a chaotic quantum phe- 
nomenon, which has no counterpart in the classical limit was given [5]. One 
of the solutions of this obstacle is based on the idea of taking the number of 
degrees of freedom to infinity [6]. For example, in [7] it was shown that a 
class of quantum dynamics of harmonic crystals becomes ergodic and mixing 
in the thermodynamic limit. Moreover, by taking ^ — > 0, classical properties 
of ergodicity and mixing are recovered. Similar results, but for the ideal gas 
quantized according to the Maxwell-Boltzmann statistics, are presented in 
[8]. 

Another way to recognize chaotic behavior in quantum systems is to inves- 
tigate the concept of entropy. Kosloff and Rice [9] introduced a generalization 
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of Kolmogorov- Sinai entropy for the quantum case, which allows to compare 
the behavior of a given system when described alternatively by classical and 
quantum mechanics. An even simpler idea was proposed by Thiele and Stone 
[10]. They suggested that the von Neumann entropy of the time-averaged 
density matrix (minus the entropy connected with the preparation of the 
initial state) can measure quantum chaos. The entropy point of view was 
taken also in a more recent paper by Slomczyhski and Zyczkowski [11]. They 
considered a pair consisting of a quantum system and a measuring appara- 
tus and proposed a new definition of entropy to measure chaotic behavior 
of such a coupled system. According to them "the approach linking chaos 
with the unpredictability of the measurement outcomes is the right one in 
the quantum case". 

On the other hand papers concerning the possible generalization of the 
notion of Lyapunov exponents have appeared. Because there is no quan- 
tum analog of the classical trajectory, the starting point has to be different. 
Perron-Frobenius operators acting on the space of densities provide a natural 
frame for the construction of the quantum counterpart of classical charac- 
teristic exponents. Such a construction has been carried out for example in 
[12,13]. 

In the present paper we will not discuss this interesting subject in general. 
Leaving aside the problem of existence and definition of an intrinsic quantum 
chaos, we adopt the point of view of [11] and restrict ourselves to the case 
of a quantum system interacting with a measuring device, or, more broadly, 
to quantum open systems. Because, as was shown by Graham [14,15], in a 
class of such systems the sensitive dependence of some expectation values on 
initial conditions remains, and the limits t oo, ^ ^ may commute for a 
dissipative dynamics, so the notion of chaotic behavior seems to make sense 
there. Moreover, in the presence of dissipation, coherence effects degrade 
and give way to an incoherent dynamics closer to the classical behavior. One 
of the attempts in this direction was the investigation of properties of the 
quantum survival probability function in open systems [16]. For the following 
master equation 

p = -i[H, p] + ^{[Hp, H] + [H, pH]) 

the different behavior in the regular and chaotic cases of the quantum sur- 
vival probability function averaged over initial conditions and Hamiltonian 
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ensembles was demonstrated. For a general discussion of quantum chaos 
with dissipation see [17]. A more general approach investigating hypercyclic- 
ity and chaos in the context of strongly continuous semigroups of bounded 
linear operators in Banach spaces was proposed in [18]. However, the def- 
inition of a chaotic semigroup given there cannot be applied to contractive 
semigroups such as Perron- Frobenius semigroups acting on the space of den- 
sities. In this case the idea of exactness of the system proved to be fruitful 
in the description of chaotic Markov semigroups associated with some differ- 
ential equations [19]. 

In this paper we generalize the notions of completely mixing and exact 
systems to the quantum level and propose a quantity, the quantum charac- 
teristic exponent Ag, which, in the classical case, corresponds to the highest 
order Lyapunov exponent measuring the speed of convergence of an exact 
classical system to statistical equilibrium. The property > selects a 
subclass of completely mixing systems which we call exponentially mixing. 
However, contrary to the classical case, exponentially mixing quantum open 
systems may not imply chaotic behavior. The relation of Xq to that one 
proposed by Majewski and Kuna in [13] is also discussed. Finally, these 
concepts are illustrated by the examples of quantum measurements based on 
event enhanced quantum theory (EEQT). 

2. Classical systems 

Classical mechanics deals with trajectories of dynamical systems whereas the 
time evolution in quantum mechanics is formulated (in Schrodinger picture) 
in terms of density matrices. They correspond to integrable positive func- 
tions in the commutative case. Therefore, to be closer to the framework of 
quantum mechanics, we formulate some classical concepts of ergodic theory 
in terms of densities. 

We start with recalling an intimate connection between the behavior of 
trajectories of a dynamical system and the evolution of its densities. For a 
general discussion of this point see [19]. We consider mainly the discrete time 
case. The reformulation of appropriate forniTilas to the case of continuous 
time systems is straightforward. Suppose (X. A. /i) is a normalized measure 
space (Ai(^) = 1) and S : X ^ X a. measure preserving transformation of 
X, that is: 

-^AeA S-\A) e A, 
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Having defined S we choose an initial point xq E X and observe its trajectory 
(xo, S{xo), S'^{xo), ...). The chaotic behavior of the system can be recognized 
by a high sensitivity of the trajectory with respect to a shght change of the 
initial state. It means that if we start with a set of initial conditions located in 
a small region then, after a large number of iterations of the transformation 5*, 
the points fill completely the whole space X. More precisely, for S"'{A) G A 

lim = 1 yAeA, ii(A) > (1) 

A transformation satisfying (1) is called exact. Notice, that invertible trans- 
formations can not be exact. A typical example of such a behavior is the 
r-adic transformation S : [0, 1) — > [0, 1) given by S{x) — rx{mod 1), 
r = 2, 3,.... 

The behavior of trajectories can be described in terms of densities i.e. 
nonnegative, integrable functions on X with integral being equal to one. 
The space of densities will be denoted by D: 

D ^ {fe L\X, ^, /.):/> a.e., = 1} 

It is well known that densities and trajectories are related in the following 
sense. 

Theorem 2.1 [19]. S is exact if and only if for all f E D 

lim IIP"/ - l||i = 0, 

n— »oo 

where P is the Perron- Frohenius operator corresponding to the transformation 
S, P : L^{X, A, jj) L^{X, A, n), given by 

j{Pf){x)dii{x) = I f{x)d^{x), ^AeA 

A S-'^{A) 

For example, in the case of r-adic transformation and for the density f{x) = 
2x, we have 



and so HP"/ — l||i = 

The concept of entropy gives a useful condition for a system to be exact. Let 

H{f) = / r){f{x))dix{x), (2) 

X 

where r]{x) = —xlogx for x > and r]{0) = 0. Then lim.„^oo = 
for all bounded f E D implies that system S is exact. Observe that is 
the maximum value of the entropy H. If the measure /j, is finite but not 
normalized, then the above theorem should be modified as follows. If 

^lim//(PV) = H^a. = H{-^yl) (3) 

then S is exact. If /u(X) = oo, then there is no constant density and we 
replace the condition (3) by the following one 

jim //(PV|P"^) = for f,geD (4) 
Here H{f\g) is the relative entropy defined by 

H{f\g) = / /(^) log J^^^^^'^ ^"PP-^ ^ ^"PP^ 

It has the following properties. H{f\g) > and H{Pf\Pg) < H{f\g). Ob- 
serve that, if there is a constant invariant density, then (4) implies (3). In 
fact, putting g — -j^^ ■ 1 we get 

H{P-f\P^g) = [-H{P-f) + i/(-^.l)]^0 

We take conditions (3) or (4) as qualitative indicators of the chaotic behavior 
of dynamical systems. In general, condition (4) gives less then (3) (see [20]): 

lim H(P''f\P''g) = ^ lim IIP"/ - P"o||i = 

n— »oo n-+oo 

Such a system, in which HP"/ — P"5'||i — > 0, n ^ oo, for any two densities, is 
called completely mixing. It is clear that every exact system is completely 
mixing. The converse is true if there is an invariant density for the Perron- 
Frobenius operator P. However, the heat equation generates a semigroup 
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which is completely mixing but has no invariant density what shows that the 
converse does not hold in general. 

To describe quantitatively the chaotic behavior one can use the Lyapunov 
characteristic exponents. The Lyapunov exponents of a given trajectory char- 
acterize the mean exponential rate of divergence of nearby trajectories. In 
m-dimensional phase space there are m (possibly nondistinct) characteristic 
exponents 

0-l{Xo) > Cr2{xo) > ... > (Tm{Xo) 

The largest one is defined by 

1 d{S^x,,S^xo) 
cAxo) = lim hmsup —log r — , 

''^'^ d{xi,xo)^o n d{xi,Xo) 

where d{x, y) is a metric on the phase space X. For example, in the r-adic 
case (t{xo) = logr and docs not depend on the initial point xq. 

Higher order exponents can be also defined [21]. The mean exponential 
growth of the phase space volume given by 

<7(-)(xo) = hm -log '^"^^'^;,^' = Ea.lxo), (5) 

where V{xq) is a m-dimensional parallelepiped with one vertex placed in xq, 
is of particular interest. In view of theorem 2.1 the growth of the phase 
space volume corresponds to decrease of any density to the constant density. 
Therefore, for completely mixing systems we define another quantity 

2 1 1 pn J pn I L 

A(/o) = liminf lim [ log — — — ] (6) 

for fo, f & Dq, Dq being some appropriately chosen dense subspace in D. 
Observe that for measure preserving fiows both a"^{xQ) and A(/o) are equal to 
zero. We now discuss the elements of the above formula more precisely. The 
minus sign refiects the property that for fixed / and /o ||-P"/ — -P"/o||i —>■ 
when n — > oo. We have put the minus sign because our objective here is to 
replace characteristic exponents, which are defined in terms of trajectories, 
by a quantity expressed in terms of densities in such a way that they would 
measure the same property of the given dynamics, and hence coincide for 
simple one dimensional systems. As will be shown in proposition 2.2, in 
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the r-adic case the Lyapunov characteristic exponent (t{xo) = logr, for any 
xo G (0, 1), indeed equals to our quantity A(l). Therefore, although this 
minus sign may seem to be artificial, it is necessary in order to describe 
the same feature of the dynamics. Moreover, we use the liminf expression 
because for some densities, even very regular ones, the action of the Perron- 
Frobenius operator may produce a constant density in a finite number of 
iterations, what would result in an infinite limit. The restriction of D to 
some subspace is a crucial point. It is known from the analysis of mixing 
systems that the decay of correlations depends on the functions chosen. As 
was shown in [22] in the case of a chaotic area preserving map on the torus 
there are examples of functions for which the correlations decay faster than 
exponentially, exponentially and only algebraically. Hence the rate of decay is 
sensitive to the choice of functions, and so the restriction to a subset of "nice" 
functions is unavoidable. Finally, observe that formula (6) is equivalent to 

1 II P" f P" f II 

A(/o) = inf lim[--log^i — ^] 

In fact, if fk is an arbitrary sequence from Dq, then the sequence = 
fo + ^^^T^ satisfies \\gk - /o||i ^ and 
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Wdk — /o||l Wfk — fo\\l 

Moreover, because we calculate n — > oo limit first, the normalization can be 
dropped, and so 

To show that formula (6) is indeed related with the Lyapunov exponents we 
calculate A(l) in the r-adic case. At first wc rcscalcd the interval [0, 1) into 
[0, 27r), replace dx by dx/2n and define S{x) = ra;(mod27r). Hence X can 
be seen as the circle S^. 

Proposition 2.2. Let Do = DnC\S^). Then = logr. 

Proof: Let / denote the Fourier transform of a function f e Do, f ^ fo, i.e. 

2tt 



f = ^le-"'^f{x)dx, keZ 
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Let P be the Perron-Frobenius operator corresponding to S. Because 



so 



r — r r 



J=0 



r-1 



j=0 

r 

and (P"/)^(A;) = /(A;r"). Hence 

IIP'^/ - < ||PV - = ll(i^V)^ - (8) 



= [2EI/(^^ 

A;=l 

On the other hand 



n\ 1 21 1/2 



11/ - > 11/ - = sup|/>)| (9) 



fc>l 



Because hm^^oo |/(^)| = 0, so ||/ — iH^oo = |/(m)| for some natural m. 
Since /' is continuous, hence integrable, so hnifc^oo ^1/(^)1 ~ 0. It imphes 
that there exists ko such that for any k > ko we have A;|/(A;)| < |/(m)|. Let 
us choose n such that > /cq- Then 

[2El/>r«)|T/^ < where A = pf^jV^ (10) 

k=l k=l 

Combining (8), (9) and (10) wc arrive at 

IIP"/ - l|Ui_ ^ A 



11/ - 1| 



Hence A(l) > logr. To finish the proof one only has to find a sequence of 
densities fk such that limfc^oo ||/fc — l||i = and 

lim lim | log ^ ,, r-^l = log r 

k^^n^oo^ n ^ ll/fc - l||i 
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It may be easily checked that the sequence fk{x) = 1 + -^{x — tt), A; e N, 
fulfills the required conditions. □ 

Completely mixing systems in which A(/) > for some / G -D will be called 
exponentially mixing. As the example of r-adic transformation shows they 
are closely related to chaotic systems. Having discussed the signatures of the 
chaotic behavior in terms of probability densities for classical systems we now 
consider the quantum case. 

3. Quantum open systems 

As we mentioned in the introduction we discuss only dissipative quantum 
systems. Dissipation in quantum theory appears by the coupling of a sys- 
tem to a reservoir. The tracing over classical variables leads to the reduced 
dynamics. Using certain approximation technics [23] it can be shown that 
the reduced dynamics is given by a Markov semigroup acting on the set 
of reduced density matrices. The Hamiltonian evolution equation is re- 
placed by the master equation. Such an evolution does not conserve energy 
and maps pure states into mixed states. For technical reason we confine 
the discussion to N-level quantum systems. Hence a dynamical semigroup 

: Mjvxiv(C) — >• MjvxAr(C) maps positive definite matrices into positive 
definite matrices and preserves the trace tr. As the composition of a Hamil- 
tonian evolution and a conditional expectation Tt is also completely positive. 

As the first indicator of the mixing property we consider the relative von 
Neumann entropy H{p\a) defined by 

H{p\a) = tr{plogp - plogcr) 

for any density matrices such that suppp C suppa. It is well known that 
H is non increasing with respect to Tt i.e. H(Ttp\Tta) < H{p\a). Following 
equation (4) we define a system to be completely mixing if 

lim H{Ttp\Tta) = (11) 

t— »oo 

holds. If a totally mixed state is T(-invariant, then (11) imphes that 
H(Ttp) Hmax- Since in the classical case it leads to the exactness of a 
system we take this property as the definition of exactness in the quantum 
case as well. Let us notice that not every dissipative system is completely 
mixing. For example, for the following master equation 

p = aipai - p 
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where p is a 2 x 2 density matrix and oi denotes the first Pauh matrix, 
the t ^ oo hmit of the relative entropy of two neighboring one-dimensional 
projectors 

_ / 1 \ _ / cos^ sin COS \ 
y y \^ sin cos sin^ j 

e (0, 7r/4), is given by \im.t-,ooH{Ttei\Tte2) = —log cos 20 ^ 0. Let us 
also point out that Hamiltonian dynamics cannot be completely mixing since 
in this case the relative entropy is constant in time. 

Because H{p\a) > — a\\l, where \\A\\i = tr\A\, so 

lim H{Ttp\Tta) = ^ lim \\Ttp - Tta\\i = (12) 

t^oo t—^oo 

Therefore, for a completely mixing system one can ask the question how fast 
the hmit \\Ttp — Tta\\i, i — > oo, tends to zero. Guided by the classical expe- 
rience (compare formula (7)) we propose the following quantity to measure 
the exponential rate of convergence 

Xg{p) = iiif lim[-ilog||Tip - T^alli] (13) 

and call it the quantum characteristic exponent. It is similar to that one 
proposed by Majewski and Kuna [13] (see also [24]). Indeed, Xq{p) — 

inf cr^p Xq{p; a) and \q{p; cr) coincides (up to the minus sign) with their for- 
mula if we notice that Tt is linear and calculate the corresponding derivative 
in a tangent direction to the space of all density matrices, i.e. for 6x such 
that trSx — 0. As in the classical case completely mixing systems with 
Xq{p) > 0, for some p, will be called exponentially mixing. Finally, let us 
point out that if a completely mixing system has a stationary density matrix 
Po, that is Tt{po) = Po for all t > 0, then po is unique and Xg does not depend 
on the choice of an initial statistical state p. In other words Xq{p) = Ag(po), 
where ^ 

Xq{po) = inf lim [--log 1 1 TtfT - po||i] 

It is also worth noting that the entropy production rate for an evolving 
reduced density matrix of an open quantum system proved to be a fruitful 
indicator of the character of dynamical behavior. As was shown in [25] the 
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classical unpredictability corresponds to the rapid entropy production on the 
Lyapunov time-scale in quantum analogs of classical systems which exhibit 
chaotic behavior. By contrast, dynamics of analogs of integrable systems 
leads to a much slower evolution towards a dynamical equilibrium of the sys- 
tem. Therefore, the rate of increase of entropy or, more generally, of decrease 
of relative entropy, can indeed distinguish between chaotic and regular quan- 
tum evolutions, in a similar way as in the classical case. 

4. Continuous quantum measurements 

EEQT is a minimal extension of quantum theory that accounts for events 
[26,27]. It postulates that a quantum measurement process is a particular 
coupling between a quantum and classical system. The time evolution of such 
a hybrid system is determined by a completely positive semigroup of linear 
operators on the space of density matrices of the total system. Moreover it 
provides the interpretation of the continuous evolution of ensembles in terms 
of a piecewise deterministic process with values in the pure state space of the 
total system. The process, after averaging, reproduces the dynamical equa- 
tion for statistical states. Such a process turned out to be unique [28] , what 
allows for deducing the algorithm generating sample histories of an individ- 
ual quantum system [29]. The sensitive dependence of sample paths of the 
associated process on the initial conditions together with fractal structure of 
the limit set of a sample path indicate chaotic properties of the corresponding 
dynamical semigroup. In this way they constitute the same counterpart for 
the evolution of density matrices as trajectories do for density functions in 
the classical case of {X) spaces. 

Let us now briefly describe the framework of EEQT. Suppose that possible 
states of the measuring apparatus C form a discrete set labeled by a = 
1, 2,...,m. The algebra of observables of C is the algebra Ac of complex 
finite sequences /„, a = 1, . . . , m. For technical reason we use Hilbert space 
language even for the description of the classical system. Thus we introduce 
an m-dimensional Hilbert space He with a fixed basis, and we realize Ac as 
the algebra of diagonal matrices F = diag(/i, . . . , fm)- Statistical states of C 
are then diagonal density matrices diag(pi, . . . ,Pm), and pure states of C are 
vectors of the fixed basis of Tic- Suppose further Q is the quantum system 
whose bounded observables form the algebra Aq of bounded operators on a 
Hilbert space Hq. We assume that Hq is finite dimensional. Pure states of Q 
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form a complex projective space CP{T-Cq) over Tiq. Statistical states of Q are 
given by non- negative density matrices p, with tr{p) = 1. The algebra At 
of observables of the total system T — Q x C is given by the tensor product 
of algebras of observables of Q and C: At = Aq® Ac- It acts on the tensor 
product Hq^Hc = ®^=iHa, where Ti^ ~ 'Hq. Thus At can be thought of as 
algebra of diagonal m x m matrices A = (aag), whose entries arc quantum 
operators: Qaa G Aq, = for a 7^ /3. Statistical states of Q x C are 
given by m X m diagonal matrices p — diag(pi, . . . , pm) whose entries are 
positive operators on TCq, with the normalization tr{p) — J2at''^{Pa) — 1- 
Duality between observables and states is provided by the expectation value 
< A >p= J2atf^i^aPa)- The coupliug of Q to C is specified by a matrix 
V — idaid), where gap are linear operators: Qa/s '■ 'Hp — > Tia- We assume 
9aa — 0. This condition expresses the simple fact: there is no dissipation 
without receiving information. The evolution equation for states is given by 
the Lindblad form 

1 

Pa = -i[H, Pa] +Yl9apPp9ai3 " Pa} (14) 

where = J20 9pa9poi and {■, ■} stands for the anticommutator. We apply 
now the above scheme to three concrete measurement processes. 
4.1 Measurement of noncommuting observables. In this example we model a 
simultaneous measurement of several noncommuting observables, like differ- 
ent spin projections [30]. In such a case we calculate the quantum exponent 
Xq for the reduced dynamics given by tracing out classical parameters, and 
show that it leads to a chaotic and fractal structure on the space of pure 
states of the quantum system. An interesting investigation of quantum chaos 
in open dynamical systems in which absorption leads to the appearance of 
a fractal set in the underlying classical phase space was presented in [31]. 
In the quantum case fractal structure appears in the Husimi functions of 
eigenstates of a non-unitary evolution operator. However, there is no appar- 
ent relation between the approach taken in Ref. [31] and our discussion of 
quantum fractal sets (on the set of pure states) arising from a continuous 
simultaneous observation of several noncommuting observables. 

The measuring apparatus consists of four yes-no polarizers correspond- 
ing to spin directions fii, i — 1, ...,4, arranged at the vortices of a regular 
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tetrahedron 



ni = (1, 0, 0), fi2 = (-^, 0, ^) 




1 _ /2 _^ 
3' V3' 3^ 



Because the quantum system is a two-state system, so the quantum algebra 
is given by 2 x 2 complex matrices. We assume it evolves according to 
Hamiltonian H = ^a^, ut > 0, where ai, a2, denote Pauli matrices 






The coupling is specified by choosing four operators Oj which correspond to 
four vectors Hi 

Qi = + Oifii ■ a) (15) 

where a e [0, 1]. Notice, that for a = 1, Oj are projection operators. The 
evolution equation is given by 

Pi = -t[H, Pi] + K^aiPjtti - Ka^pi (16) 

j 

what implies the following master equation for the reduced density matrix 
P ^T,iPi 

p — —i[H, p] + K^^ttipUi — Ktt^p (17) 

i 

where > is the coupling constant and = I]i • Eq.(17) implies that a 
projection valued measure corresponding to a sharp measurement has been 
replaced by a positive operator valued measure. Clearly, the totally mixed 
state J/2 is stationary. At first we show that the dynamics given by eq.(17) 
leads to an exponentially mixing system. 

To solve (17) we assume the initial state is a pure one i.e. p(0) = |(/ + 
mo • (?), mo = (ci, C2, C3) with cf + + c| = 1. Because the dynamical 
semigroup preserves positivity and trace so a general solution is of the form 
p(t) = |(/ + m{t) ■ (?) with ||m(t)|| < 1. By direct calculations we obtain 
the following system of differential equations 

mi = — a;m2 — -na mi 
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4 2 

1712 — oumi — -na m2 

O 

ma = --^f^oi ma 

and so 

4 2-/- 

= (— C2sina;t + ci cosa;t)e~3'^" 

4 2 + 

?7i2(i) = (cisma;i + C2 cos ut)e~ 3 * 

is the solution with the initial condition m(0) = mo- Because 

\\m - = = ||m(t)|| = e"!-^* 

and the convergence above does not depend on the choice of mo, so 

mi 



A,(^/) = inf lim[-^log||p(t) - = ^Ka' (18) 



We describe now a sample path of the process associated with the dynamical 
semigroup determined by eq.(17). Assume that at time t — the quantum 
system is in the state r(0) G S*^ ( we identify here the space of pure states of 
the quantum system with a two-dimensional sphere with radius 1). Under 
the time evolution it evolves to the state f{t) which is given by the rotation 
of r (0) with respect to z-axis. Then, at time ti a jump occurs. The time rate 
of jumps is governed by a homogeneous Poisson process with rate k. When 
jumping f{t) moves to 



(1 - a'^)f{t) + 2a{l + af{t) ■ ni)ni 



with probability 



Pi{r{t)) 



1 + q;2 + 2ar{t) ■ Ui 
1 + + 2ar{t) ■ Hi 



4(1 + q;2) 

And the process starts again. The iterations lead to a self-similar structure 
with sensitive dependence on the initial state. Fig. 1, 2 and 3 depict sample 
paths of the quantum particle for different values of a in the case when 
H = and k, = 1. Fig. 4 a-c present a = 0.7 case together with its zooms 
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x2 and x4. 

Numerical simulations performed for H = and k = 1 [32] show that when 
a increases from 0.75 to 0.95, then the Hausdorff dimension of the limit set 
decreases from 1.44 to 0.49. It is worth noting that at the same time the 
quantum characteristic exponent increases as «^ (see formula (18)). This 
suggests an intimate relation between and the Hausdorff dimension of the 
limit set. In the classical case such a relation was established in [33,34]. 
4.2 Quantum Zeno effect. A phenomenon of keeping a quantum state from 
evolving by performing a sequence of frequent measurements was discussed 
by Misra and Sudarshan [35] and named quantum Zeno effect (paradox). 
Roughly speaking it means that the constant observation freeze the quantum 
system in its initial state. Such a situation was described, for example, in [36] 
for the case of the neutron placed in a static magnetic field and interacting 
with a device which selects one component of its spin. By this example we 
demonstrate a non-trivial dependence of the quantum characteristic exponent 
\q on the value of the coupling constant. Because the quantum system is 
a two-state system, so the quantum algebra is just the algebra of 2 x 2 
complex matrices, while the classical system being a yes-no device consists 
of two distinct points. The evolution equation for the reduced density matrix 
P e M2X2, P > 0, trp = 1, is given by 

P = -i[H, p] + K{epe - ^{e, p}), (19) 

where Hamiltonian H = and the matrix e = |(/ + <Ji) is a projector 
onto the first eigenvector of the Pauli matrix ai. The coupling constant 
K measures the frequency of "checking" whether the quantum system is in 
the eigenstate of ci or not. Again, as in the previous example, the totally 
mixed state 1/2 is stationary. Now eq.(19) leads to the following system of 
differential equations 

Xi — —ijJXi 
±2 = OJXx - -X2 (20) 

where p(t) = |(/ -|- x ■ a). For a given initial condition p(0) it is a routine 
exercise to solve (20), and hence to obtain the following result. 
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Theorem 4.1. Let n = Aau, a > 0. Then 

if a e [0, 1] 
if a>l 

It follows that the degree of mixing (measured in terms of A^) increases from 
zero to the maximum value uj when the frequency of measurements increases 
to the critical value 4cj. The further increase of number of observations results 
in the decrease of A„ to like -. Hence we can conclude that the interaction 
with the measuring apparatus makes the evolution of the quantum system 
exact (and so mixing), but constant observation suppresses this property. 
4.3 Two-level atom driven by a laser. Let us consider the fluorescent photons 
emitted by a single, two- level atom that is coherently driven by an external 
electromagnetic field. It is known that the quantum system evolves from the 
ground state in a dissipative way. When a photoelectric count is recorded 
by a photoelectric detector (we assume the detector efficiency to be equal to 
one), the atom returns to the ground state with the emission of one photon. 
The master equation for such a system is given by the following formula 

p = -i[H, p] + ^ApA* - p} (21) 

where 

n 

H — (7i and A 

2 

Vt is Rabi frequency and 7 > is the relaxation rate. Because the unique 
stationary state po is given by po = |(-^ + ^2^2 + 't-sCs), where 

2O7 -72 




4a;2 + ^2 ^ 4^2 + ^2 

so the system is not exact. However, it is exponentially mixing. Solving 
eq.(21) and using formula (13) we obtain that A^(po) = 7/2. Hence, al- 
though Ag(po) > it would be difficult to consider this system chaotic. The 
property of being exponentially mixing is a new feature in quantum open 
systems. 
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